Non-Fermi-Liquid Transport Phenomena and Superconductivity Driven by Orbital 
Fluctuations in Iron Pnictides: Analysis by Fluctuation- Exchange Approximation 
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We study the five-orbital Hubbard model including the charge quadrupole interaction for iron 
pnictides. Using the fluctuation-exchange approximation, orbital fluctuations evolve inversely pro- 
portional to the temperature, and therefore the resistivity shows linear or convex T-dependence for a 
wide range of temperatures. We also analyze the Eliashberg gap equation, and show that an s-wave 
superconducting state without sign reversal (s++-wave state) emerges when the orbital fluctuations 
dominate over the spin fluctuations. When both fluctuations are comparable, their competition 
gives rise to a nodal s-wave state. The present study offers us a unified explanation for both the 
normal and superconducting states. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 
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I. INTRODUCTION 

The many-body electronic states and the pairing mech- 
anism in iron pnictides have been significant open prob- 
lems. By taking account of the Coulomb interaction and 
the nesting of the Fermi surfaces (FSs) in Fig. 1(a), a 
fully-gapped sign-reversing s-wave state (s±-wave state) 
had been proposed ^"^. Experimentally, both Tc and an- 
tiferro (AF) spin correlation increases as x decreases in 
BaFe2(Asi_xPx)2^- In contrast, Tc in LaFcAsOi-^^Fa; 
at X = 0.14 increases from 23 K to 43 K by apply- 
ing the pressure, whereas AF spin correlation is almost 
unchanged^. Thus, the relationship between Tc and 
strength of the spin fluctuation seems to depend on com- 
pounds. 

On the other hand, an orbital-fluctuation-mediated s- 
wave state without sign reversal (s++-wave state) had 
been proposed based on the five-orbital Hubbard model 
including the charge quadrupole interaction®^^^. The 
charge quadrupole interaction is induced by the vertex 
correction (VC)^ due to the Coulomb interaction and the 
electron-phonon (e-ph) interaction due to Fe-ion Einstein 
oscillations. Within the random-phasc-approximation 
(RPA), it was found that d-orbital fluctuation is in- 
duced by small e-ph interaction. Especially, the em- 
pirical relationship between Tc and the As-Fe-As bond 
angle (Lee plot)^^ has been naturally explained. Re- 
cently, theoretically predicated orbital fluctuations^' ^'^ 
have been detected via the substantial softening of the 
shear modulus^^. The softening of the shear modulus 
and the structure transition have been explained by the 
two-orbiton mechanism based on the orbital fluctuation 
theory^^. The s++-wave state is consistent with the ro- 
bustness of Tc against randomness ^^"-'^^ as well as the 
"resonance-like" peak structure in the neutron inelastic 
scattering^^. 

However, spin/orbital fluctuations obtained by the 
RPA are reduced by the self-energy correction. There- 
fore, in order to confirm the orbital fluctuation sce- 
nario, it is desired to analyze the many-body electronic 



states beyond the RPA. For this purpose, the fluctuation- 
exchange (FLEX) approximation^^ would be appropri- 
ate, in which the absence of spin/orbital order in 2D 
systems, known as the Mermin- Wagner theorem, is rig- 
orously satisfied^°. 

In this paper, we analyze the five-orbital Hubbard 
model including the charge quadrupole interaction for 
iron pnictides using the FLEX approximation.^^ In the 
normal state, large orbital fluctuations induce highly 
anisotropic quasiparticle lifetime on the FSs as well as the 
T-linear or T-convex resistivity ^22-24 g^j-^j ^j^g large neg- 
ative thermo-electric power S. The large orbital fluctua- 
tions also introduce the S-|.+-wave superconducting (SC) 
state for a wide range of parameters, and the competition 
between orbital and spin fluctuations lead to the nodal 
s-wave state. We propose that the orbital fluctuation is 
the origin of both the s++-wave SC state and the non- 
Fermi-liquid transport phenomena in the normal state. 



II. FORMULATION 

In this paper, we set the x and y axes parallel to the 
nearest Fe-Fe bonds and the orbital z^, xz, yz, xy and 
x^ — y^ orbitals are denoted as 1, 2, 3, 4 and 5, respec- 
tively. 

We employ the five-orbital Hubbard model^ including 
the quadrupole-quadrupole [electron-electron (el-el)] in- 
teraction induced by the VC due to the Coulomb interac- 
tion and e-ph interaction due to Fe-ion Einstein optical 
modes. The quadrupole-quadrupole interaction is given 



as 
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V{Ur, 



~g(w„) Y. [OIM. + 01 A + OlyOly) ,(1) 



where Op (F = xz,yz,xy) is the charge quadrupole 



operator and g{(jJn) = gujf)/{uj'^ -f Wp) is proportional 
to the phonon Green function; g — g{Q) is the effec- 
tive el-el interaction for a;„ — 0, and ojd is the cutoff 



frequency^. For example, we show non-zero Vu 



for 



1,1', m,m' = 2,3,4 in Fig. 1(b). Other than Fig. 1(b), 
V has many non-zero off-diagonal elements as explained 
in Ref."'^°'^^, since the Fe-ion oscillation (non-Aig mode) 
induces various inter-orbital transitions. 

In the FLEX approximation^^, the 5x5 self-energy 
matrix E in the orbital representation is given by 

^l^lM ^^Y.T. ^hi.MiM)Gi.iAk - q), (2) 
q hii 

where li represents the orbital, A'^ is the number of k 
meshes, and we denote k = [k, tn) with fermion Mat- 
subara frequency e„ = (2n -f- l)7rT, and q = (q,ujn) with 
boson Matsubara frequency lo„ = 2mTT. 

G is the 5x5 Green function matrix in the orbital basis, 
and V^ is the 25 x 25 interaction term for the self-energy 
given as^^ 

V'^iq) = \t^r{q)^' + \^'r{q)^' 

- \{t^- f ^)x'"(«)(f ^ - f^) + ^f ^- + ^f^(3) 
where the irreducible susceptibility is given by 

^"i..hiM) = -^Y.^hhik + q)GuiM. (4) 

k 

and the spin (orbital) susceptibility is obtained as 



^^(0 = 



l-r"(c)yi 



X" 



(5) 



Here, F" = 5 [F^ = -C - 2F(a;„)] is the irreducible 
vertex for the spin [charge] channel; S and C repre- 
sent the Coulomb interaction in the multiorbital model 
introduced in Refs.^'^'^°'^^; Their matrix elements con- 
sist of the intra-orbital Coulomb U , the intcr-orbital 
Coulomb [/', Hund's coupling J and the pair hopping 
J' . We assume that 3 = 3' and U = U' -\- 23. Since 
the Fe-ion oscillation induces various inter-orbital tran- 
sitions, the substantial orbital fluctuations appear at 
low frequencies. On the other hand, the charge sus- 
ceptibility x^iq) = EimXH,mm(9) ^^ "ot enhanced due 
to the cancellation^'^". In the present study, we drop 
ladder-type diagrams by V{uin), which is justified when 
wd ^ £'f^'^°. For the same reason, V{ujn) is absent in 

In the FLEX approximation, we obtain G and E self- 
consistently using the Dyson equation G~^ = (G^)^^ — E. 
In multiband systems, the FSs are modified from the orig- 
inal FSs due to the self-energy correction. To escape from 
this difficulty, we subtract the constant term [E(fc, -)~iO) + 
E(fc, — i0)]/2 from the original self-energy, corresponding 
to the elimination of double-counting terms between LDA 
and FLEX26. Hereafter, we fix 3/U = 1/6, wq = 0.02eV, 
and the electron filling n = 6.1 except for Fig. 5. Because 



of the smallness of the FSs in Fig. 1, fine k meshes are re- 
quired for a quantitative study. We take N = 128 x 128 
k meshes which is four times that used in Ref.^^, and 
1024 Matsubara frequencies. Then, we obtain reliable 
numerical results for T > O.OleV. 
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FIG. 1: (color online) (a) FSs in the unfolded zone. The dot- 
ted circles represent the cold-spot given by the orbital fluc- 
tuation theory. The cold-spot is composed of xz/yz-orhitals. 
(b) Phonon-mediated el-el interaction (V) for 2,3,4 orbitals. 
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FIG. 2: (color online) (a) fe-dependence of ^"'^'^' induced by 
the spin (orbital) fluctuations on each FS. Note that 7" de- 
creases with g due to the suppression of x" by 7^^. (b) T- 
dependence of Si,(c) = (1 ~ Q;s(c))~^- (c),(d) T-dependence 
of p. p = 1 corresponds to {hac)/e^ ~ SOO/icmfi when the 
interlayer distance is Oc = 0.6nm. 



III. RESULT 

A. Normal state 

We begin with the electronic property in the normal 
state. Hereafter, the unit of energy is eV. First, we dis- 
cuss the quasiparticle damping rate 7^ on each FS, which 
is given by the imaginary part of the self-energy in the 



band-diagonal representation. In Fig. 2(a), 7^ repre- 
sents the damping due to spin (orbital) fluctuations for 
T = 0.015 and U = 1.8, which is given by substituting 
V^ = |f"x'r" + |f" (^f ^x'^r^ + h^") in Eq. (2). The 
horizontal axis is the azimuth angle for the fe point with 
the origin at the r(M) point for FS1,2 (FS4). The re- 
lationship 7fc ~ 7fc -f 7fe is satisfied since the third term 
in Eq. (3) is very small. We will see below that the 
value [/ = 1.8 can reproduce moderate AF spin fluctua- 
tions observed in e-doped compounds, and it is consistent 
with C/ ~ 2 reported by x-ray absorption spectroscopy 
(XAS)27. 

In Fig. 2(a), the relation 7'' ':^ 7^ holds for 5 = 0, and 
the momentum dependence of 7^ on each FS is small al- 
though the AF spin correlation is well developed. The 
value of 7^ increases with g, and 7"^ ^^ 7" at 5 = 0.26. 
In Fig. 2(a), 7^^ ^ 7^* for g — 0.3; the corresponding 
dimensionless coupling is just A = gN{0) ^ 0.2. Then, 
7£ on FS4 (e-pocket) is anisotropic due to the orbital de- 
pendence of xii',mm', and it takes the minimum value at 
~ 0, where the FS is composed of 2, 3-orbitals^. This 
"cold-spot" is important for the transport phenomena. 
Since the cold spot is on the e-pocket, the Hall coefficient 
Rn will be negative, consistent with experiments^^'^^'^^. 
In the case of high-Tc cuprates, various non-Fermi-liquid 
transport phenomena (e.g., violation of Kohler's rule) 
originate from the cold/hot spot structure as well as the 
backflow (=current vertex correction) due to the spin 
fluctuations'^^. By analogy, the appearance of the cold 
spot in Fig. 2(a) indicates that the orbital fluctuations 
are the origin of striking non- Fermi-liquid transport phe- 
nomena in iron pnictides^^'^'^^^^. 

In Fig. 2(b), we show how the orbital and spin fluc- 
tuations develop as T decreases: In the FLEX, the 
spin (orbital) susceptibility is enhanced by the spin (or- 
bital) Stoner enhancement factor Ss(c) = (1 ^ <^s(c))~^j 
where agtc) is the maximum of the largest eigenvalue of 
ps(c)^irr^q^Q^ with respect to q. Then, as^c = 1 corre- 
sponds to the spin/orbital order, although it is prohibited 
in 2D systems by the Mermin- Wagner theorem^°. In the 
case of C/ = 1.8 and 5 = 0, large Sg (^ 10) is produced 
at q w Q = (7r,0) (i.e., x'(Q,0) ex Ss)- Ss gradually 
increases as T drops, which is a typical critical behavior 
near the AF magnetic quantum-critical-point (QCP)^°. 
When g > 0, x^(9: 0) is enhanced at q = and q = Q 
almost equivalcntly^". At g — 0.3, large Sc (^ 10) is pro- 
duced at <7 pa Q or 0, and it increases approximately pro- 
portional to T~^. Thus, it is confirmed that both ferro- 
and AF-orbital fluctuations show critical evolutions near 
the orbital QCP. 

Next, we discuss the resistivity p due to the orbital 
and spin fluctuations. By neglecting the backflow, the 
conductivity is obtained by 



N '^ 

k.a 



duj 



"a.k 



GkA^ + iO)\ ,(6) 



index, /(w) is the Fermi distribution function, v^ ,^ is the 
velocity of band a, and Gk^a{i^+iO) is the retarded Green 
function for band a in the FLEX approximation. Figure 
2(c) and (d) show the obtained resistivity p = 1/axx for 
U — 1.2 and 1.8: In case oi U = 1.2, p shows a conven- 
tional sublinear (concave) T-dependence at 5 = 0. p in- 
creases with g due to the orbital fluctuations, and almost 
T-linear resistivity is realized a.t g — 0.22. At g — 0.25, 
p shows a superlinear (convex) T-dependence. In case 
of [/ = 1.8, p is linear-in-T at g = 0, while it shows 
a clear superlinear T-dependence at g — 0.3. Thus, we 
stress that non- Fermi-liquid resistivity is not a direct evi- 
dence for the spin fluctuations. In LnFeAsO compounds, 
Tc increases as the AS4 tetrahedron is close to a regular 
one, and the T-dependence of p changes from concave 
to convex^^. Since g is maximum when the As4 tetrahe- 
dron is regular^*^, this experimental correlation between 
Tc and p(T) is understood in terms of the orbital fluc- 
tuation scenario. We note that non-Fermi-liquid-like fre- 
quency dependence of ImS(w) was recently discussed in 
Ref.3i. 
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where e(< 0) is the charge of an electron, a is the band 



FIG. 3: (Color online) T-dependence of 5* for (a) U = 1.2 and 

(b) (7=1.8 

Here, we study the non-Fermi-liquid-like behavior of 
the thermo-electric power S induced by the orbital and 
spin fluctuations. According to experimental results"^^"^^, 
S in e-doped systems is negative below the room tem- 
perature, and \S\ develops inversely proportional to 
the temperature till T* ^ lOOK. In optimum doped 
Ba(Fei_2;Co2;)2As2 peak value of \S\ ^ 50/iV/K and 



T* ^ 130K are observed. '^^ Similar behaviors are also 
observed in high-Tc cuprates with the strong AF spin 
fluctuation. Thus, such remarkable non-Fermi-liquid be- 
haviors in iron pnictides are expected to be realized by 
the orbital and spin fluctuations. 

Since the effect of backflow is small for S,^"^ we calcu- 
late S by neglecting the backflow as follows 



5=- 



Ac, a 



duj f df{uj) 
duj 



^ \v% kGk,a{^ + iO)\ 



The obtained T-dcpendence of S is shown in Fig. 3. In 
the case oi U = 1.2, Fermi-liquid-like behavior {S oc T) 
is obtained for 5 = 0, where both the spin and orbital 
fluctuations are weak. For g = 0.25, where the orbital 
fluctuation is strong, \S\ becomes larger than that for g = 
0, and the deviation from the Fermi-liquid-like behavior 
is realized. In the case oiU = 1.8, where spin fluctuation 
is strong, non-Fermi-liquid-like behavior is obtained as 
shown in Fig. 3(b). For g = 0, the value of \S\ is small 
and almost independent of T. On the other hand, 15*1 is 
drastically enhanced, and shows the peak at T* ~ 150K 
for g = 0.3, where both the spin and orbital fluctuation 
are strong. The obtained result for [/ = 1.8 and g = 0.3 
is consistent with experiments. '^^"'^^ In the following, we 
explain why the absolute value of S becomes large for 
large value of g. Since Va^k ~ ^/Na{uj) at w = e^, where 
Na and e^ are the density of state and the dispersion on 
band a, respectively, S is rewritten as 
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E 

a 

E^ 
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Zr.UJ 



dfiio) 



d_ / 1 

duj V^a(w)7o(w)/<^=o 



duj 



(8) 



where 



= 1- 






and 7a are the renormal- 



ization factor and the quasiparticle damping on band a, 
respectively. When the orbital fluctuations are weak, S 
takes small and negative value because -§^Na > is sat- 
isfied on the e-pocket (cold spots) . In the case of strong 
orbital fluctuation with large value of g, S" is still neg- 
ative while the absolute value is much enhanced, due 
to the large value of ^7q > at the cold spots on 
e-pocket shown in Fig. 2(a). Thus, the orbital fluctu- 
ation plays an important role in enhancing the absolute 
value of iS* and reproducing the experimental results. We 
stress that result for U = 1.8 and g = 0.3 well repro- 
duce the experimental behaviors of S in optimum doped 
Ba(Fei„,Co,)2As2.36 

In the present FLEX approximation, almost isotropic 
damping is obtained for U = 1.8 and 5 = 0, as shown in 
Fig. 2(a). In contrast, Kemper et al.^^ reported a clear 
hot/cold spot structure due to the spin fluctuation using 
the RPA self-energy (self-inconsistent FLEX approxima- 
tion). This difference would come from the presence (ab- 
sence) of self-consistency in the former (latter) calcula- 
tion. Kemper et al.^^ also reported interesting doping de- 



pendence of the sign of Hall coefficient i?H- However, cur- 
rent vertex corrections would be necessary to reproduce 
the magnitude and T-dependence of i?H appropriately.^^ 
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FIG. 4: (Color online) (a) U-g phase diagram given by solving 
the linearized Eliashberg equation at T = 0.015. Nodal s- 
wave gap state is obtained in the shaded area for nimp = 0.02 
and 0.05, and solid lines (dotted lines) represent the boundary 
between fully-gapped s+_-wave (s++-wave) state. Dashed- 
dotted line denotes Oc — 0.98. (b) s++-wave gap {Xe = 0.59) 
for C = and g = 0.24. (c) s+_-wave gap {Xe = 0.49) for 
U — 1.8 and g = and. (d) Nodal s-wave gap {Xe ~ 0.28) 
for U = 1.2 and g = 0.15. 



B. SC state 

Next, we discuss the SC state. In the presence of dilute 
impurities (nimp ^1), the linearized Eliashberg equation 
in the orbital basis is^ 



AEA„,(fc) =- ^ E W^i™i,w(fc-fc')G:„,™.(fc') 

X A„,™3 (fc')G;,,„,3 (-fc') + SJ:t, (e„),(9) 

where Aw (fc) is the gap function and Ae is the eigen- 
value that reaches unity at T = Tc. (5S° represents the 
impurity-induced gap function. (G')^^ = i^)^^ ^ ST,^, 
where G is the Green function given by Eq. (2), and (5S" 
is the impurity-induced normal self-energy. The pairing 
interaction W in Eq. (9) is 



w{q) = Ir^nqW' 



1 



r^r(9)^ 



1, 



1 



r", (10) 



where x*''^ is given by the FLEX approximation for 
^imp = 0, because of the fact that the fully self-consistent 
FLEX with impurity-induced self-energy leads to un- 
physical reduction in x^ , unless vertex correction is taken 
into account^^. The first (second) term in Eq. (10) works 
to set Afsi,2 • Afs3,4 < (> 0). 

In the T-matrix approximation, JS"'" is given as 






'T-impT'y(en), (H) 

,^T.a{en)fhn{£n)Tjm{-en), (12) 



^inip 



Im 



where Tij{en) = /(I — /g(e„)) ^ is the T-matrix in the 
normal state^^; g(e,i) = j^ J2k Gki^n) is the local normal 
Green function, and / is the local impurity potential that 
is diagonal in the orbital basis. We put 1=1 hereafter. 
In Eq._ (12), /,,(e„) = ^ Efc.;™ G,,(fc)Az„,(/c)G,>„(-fc) 
is the linearized local anomalous Green function. 

In Fig. 4(a), we show the U-g phase diagram obtained 
by the FLEX approximation. The dashed-dotted line 
represents the condition ac = 0.98 at T = 0.015, corre- 
sponding to g = 0.25 ^ 0.3. (In the RPA, the same con- 
dition is satisfied for g = 0.21 ~ 0.23.^°) Therefore, sub- 
stantial orbital fiuctuations emerge for A = gN{0) < 0.2 
even if the self-energy correction is taken into account. 
On the other hand, a, = 0.95 (0.92) for U = 1.8 
and g = (0.3) in the FLEX approximation, although 
t/cr = 1-25 for as = 1 in the RPA. Thus, the renormal- 
ization in as is much larger than that in ac, because of 
the difference in the coefficients (in factor 3) between the 
first and the second terms in Eq. (3). 

Next, we solve Eq. (9) with high accuracy using the 
Lanczos method at T = 0.015. Then, the s+^-wave gap 
function is obtained around the line ac = 0.98; Figure 
4(b) shows the s++-wave gap for g = 0.24 and U = 
{Xe = 0.59). On the other hand, s±-wave gap is ob- 
tained when g is sufficiently small; Figure 4(c) shows 
the s±-wave gap for [/ = 1.8 and .g = {Xe = 0.49). 
When riinip = 0, the gap function changes from (b) to 
(c) discontinuously on the phase boundary in Fig. 4(a), 
as found in Ref.^*^. When riimp > 0.02, however, the gap 
function changes smoothly during the crossover. Then, 
line-nodes inevitably appear on FS3,4 in the shaded area 
in Fig. 4(a); Figure 4(d) shows the nodal s-wave gap for 
U = 1.2, g = 0.15 and n^^p = 0.02 (A^ = 0.28). Thus, 
both regions for s-|-+-wave and nodal s-wave states are 
extended by a small amount of impurities, although Xe 
for the latter state is reduced by impurities. A nodal 
s-wave solution at ni,„p = with larger Xe may be ob- 
tained by considering a 3D nodal-line structure in a 3D 
tight-binding model***^. 

Here, we discuss that line nodes originate from the 
competition between the orbital and spin fluctuations: 
The electrons at 6 ^ {tt/2) on FS4 are composed 
of orbital 2,3 (4). Since the orbital 4 is absent in 
FS1,2, the nesting-driven spin correlation between the 
orbital 2,3 on FS1,2 and the orbital 4 on FS3,4 is weak. 
(That is, X24,42(Q) "^ X22,22(Q)-) On the other hand. 



both X24 42(9) ^'^d X22 22(9) dcvclop wcll^"^^. There- 
fore, when orbital and spin fluctuations are comparable, 
Afsi,2 ■ Afs4 is negative (positive) at ~ {n/2) due to 
the orbital-dependences of the spin and orbital suscepti- 
bilities. 
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FIG. 5: (Color online) (a) n dependence of Ae for s++- and 
s+_-wave states at T = 0.02 and nimp = 0. (b) g dependence 



of Ae for s++-wave state at T = 0.015 and nin 



0. 



In Fig. 5(a), we show the filling dependence of Ae for 
the s++-wave state {U = 1.2, g = 0.25), and that for the 
s+_-wave state {U = 1.8, g = 0). We note that FS1,2 
disappear for n > 6.3. The value of Ae for the s++-wave 
state decreases monotonically with n, while Ae for the 

s^i -wave state is rather insensitive to n, maybe because 

the temperature, T = 0.02, is rather high. Figure 5(b) 
shows that Ae for the s_|__|_-wave state {U = 1.2, n = 6.1) 
increases with g. 



IV. CONCLUSION 

We performed the FLEX approximation in the multi- 
orbital Hubbard model including the charge quadrupole 
interaction for iron pnictides. It was confirmed that the 
orbital-fluctuation-mediated s+4.-wave state is realized 
by small e-ph interaction g. As increasing the value of g, 
both the Tc of the s+^^-wave state and the resistivity p 
are increased, and the latter changes from T-concave to 
T-convex. This correlation between T^. and p is consis- 
tent with experiment^^. Moreover, the obtained thermo- 
electric power S" is a large negative value due to cold 
spots on the e-pocket, when the orbital fluctuation is 
dominant. The large negative value of S is consistent 
with experiments'^^"'^^. We note that the region of s++- 
wave or nodal s-wave states is enlarged in the presence 
of a small amount of impurities. Thus, the present or- 
bital fluctuation scenario presents a unifled explanation 
for both normal and SC electronic states. 
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